function plot_syn_diff(connectfile, synfile0,synfile1,neuron_num);

%connectfile='connect.dat';
[a b c d]=textread(connectfile,'%f%f%f%f%*[^\n]','delimiter','\t');
connect=[a b d]; % x y coord of area & d is cell type index.

X7a=load(synfile0);
X7b=load(synfile1);

%neuron_num=43511+16

for kk=1:length(neuron_num)
    figure(kk)
    ind=find(X7a(:,1)==neuron_num(kk)+4800+96);
    length(ind)
    syn_diff=X7b(ind,3)-X7a(ind,3);
    hold off
    for ii=1:length(ind)
        % xyz=[connect(X7_19350(ind(ii),2)+1,[1,2,8]);connect(X7_19350(ind(ii),1)+1,[1,2,8])];
        xyz0=connect(X7b(ind(ii),2)+1,:);
        if syn_diff(ii)>0.2
            plot3(xyz0(:,1),xyz0(:,2),xyz0(:,3)+1,'r.');
        elseif syn_diff(ii)<-0.2
             plot3(xyz0(:,1),xyz0(:,2),xyz0(:,3)+1,'b.');
        end;

        hold on
    end;
    xyz1=connect(X7a(ind(1),1)+1,:);
    plot3(xyz1(:,1),xyz1(:,2),xyz1(:,3)+1,'k*');
    
    title(['Neurons connected to cell #: * ', num2str(neuron_num(kk))]);
    pause;
    clear ind
end;





  